Mechanical cat states in graphene resonators 
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We study the quantum dynamics of a symmetric nanomechanical graphene resonator with degen- 
erate flexural modes. Applying voltage pulses to two back gates, flexural vibrations of the membrane 
can be selectively actuated and manipulated. For graphene, nonlinear response becomes important 
already for amplitudes comparable to the magnitude of zero point fluctuations. We show, using 
analytical and numerical methods, that this allows for creation of cat-like superpositions of coherent 
states as well as superpositions of coherent cat-like non-product states. 



Coherent superposition of states are characteristic 
traits of quantum mechanics. These phenomena have 
already been realized in many-particle contexts such as 
trapped ultra-cold atoms [1], superconductors [H and 
photonic systems [3]. A current challenge is to observe 
these effects for collective degrees of freedom in a macro- 
scopic context in, e.g., mechanical resonators Q. 

Recent advances in cooling mechanical resonators and 
sensitive displacement detection have allowed reaching 
the motional ground state and observing zero point fluc- 
tuations of center of mass [5, 6]. Active manipulation and 
characterization of the quantum state of these systems, 
as already achieved with photons [7], seem to be within 
reach. For a mechanical system, a desirable state to gen- 
erate is a 'cat' state. This is a coherent superposition 
of two minimum uncertainty wave packets separated by 
more than their individual quantum fluctuations. 

For the harmonic oscillator a minimum uncertainty 
wave packet is a coherent state \a) = exp[aa''" — a*a]JO) 



generated by displacing the oscillator ground state 
As shown by Yurke and Stoler for a nonlinear os- 
cillator H = fiuoQU + fiVtii? an initial coherent state \a) 
will after a time t = 7r/{2Q) evolve into the cat state 
{l/^/2) [e~*^/^ I a) + e*^/^ with the maximum spa- 

tial separation A = 2\a\. 

Nanoelectromechanical resonators are typically intrin- 
sically nonlinear The amplitudes needed to observe 
nonlinear effects are often orders of magnitude larger 
than the quantum zero point fluctuations xq = \/h/muJo. 
A cat state obtained due to this nonlinearity would have 
a separation A ^ xq. As the decoherence rate scales 
as (A/xo)^ [III, im this has, until now, been unfeasible. 
Instead, coupling to auxiliary quantum systems has been 
proposed to engineer the nonlinearity |13L Il4|. 

We show how the intrinsic nonlinearity in a graphene 
membrane resonator can be used to prepare cat states by 
applying voltage pulses to local backgates. The reason 
for using graphene is the ultra-low mass of the graphene 
sheet which leads to a large xq, and an onset of non- 
linear response at small amplitudes 15, [l6|. This im- 
plies that cat states with moderate ratios (A/xq) can 
be constructed without the need for engineered auxiliary 
quantum systems or feedback loops. Another feature of 



FIG. 1: (Color online) Cross section of a square graphene 
membrane resonator. A fully clamped graphene membrane 
with side L is suspended a distance uo above the substrate. 
Below, covering two adjacent quadrants beneath the mem- 
brane, are local backgates with time dependent voltage biases 
Vi,2{t). By applying pulses to the local gates, cat-like super- 
positions of flexural mode states, as well as superpositions of 
coherent cat-like non-product states, can be generated. 



two-dimensional mechanical resonators is that they can 
be designed to have degenerate flexural modes. Coupling 
between these modes can be controlled by external per- 
turbations such as gate electrodes, which can be utilized 
for state manipulation. 

For concreteness we consider a square graphene mem- 
brane with mass density po = 7.6 • 10~^ kg/m^ and side 
length L. The sheet is suspended in the x?/-plane at a 
distance uq above two local backgates (see Fig.[T]). They 
cover adjacent quadrants below the membrane and are 
have voltages Vi{t) and V2{t). The Hamiltonian den- 
sity of the system can be divided into two parts: = 
^(x) + .^^(x, Vi, where x = (x^y). Here ^(x) 
gives the intrinsic mechanics of the membrane while the 
coupling to the gates is described in .^^(x, Fi, V2). To a 



first approximation one finds [15|, [2 
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where ii(x, t) is the out of plane displacement and ttq 
its conjugate momentum density. The built-in tension is 
To and the stretching-induced tension is determined by 
Ti = (/i + A/2) where the Lame-parameters of graphene 
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are /i ?^ 3A 9 eV/A^ We model ^ in the local 
approximation as = eoT^(x, t)^/['Uo + '^(x)], where 
eo = 8.854 • 10~^^ F/m and V^(x, t) is the local gate po- 
tential. Expanding u and ttq in mode functions as 

u = ^(7kW^k(x), = ^pkW^k(x), 

k k 

using (/9k = 2 sin(/ca;x) sin(/c2yy) and k = {kx^ky) = 
(n, m)7r/L for n, m = 1,2,..., gives H = Hq + Hq with 

^0 = E ^ (f'k + A^-^c^^g^) + ^ ^ Fkk' (gk^k' 

k kk' 

Hg=Ko^Y1 ^k^k + ^ ^kk'^k^k' • 
k kk' 

Here M = poL'^ and cjk = (To/po)^|k|. The coefficients 
F^^' and i^kk' ^^^^ t)e discussed below. 

We restrict attention to the two lowest degenerate 
modes with ki = (tt/L, 27r/L), and k2 = (27r/L,7r/L) 
and their frequencies uji^2 = {tt / L) ^/ 5To / po . We label 
these modes 1 and 2. Quantizing, by imposing the com- 
mutation relation [gk,Pk'] = ^^^k,ks yields 

+/s:i(2i + K2q2 + Ki2qiq2 + ^i?? + S2ql (2) 
Introducing k = cqL'^/uqTT^ the coefficients in ^ are 



D = 



leiTr'^Ti 
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Ku = -^{V,^ - Vi), ^1,2 = -mr^VySuo (3) 

For time dependent Vi^2{t) the Hamiltonian (|2]) describes 
the excitation and evolution of two (nearly) degenerate 
interacting flexural modes. Also other modes, not in- 
cluded in (j2j), will be excited by the gates. This two- 
mode approximation is valid for weak intermode inter- 
action and when the other modes are off resonance with 
the modes 1 and 2. 

For cat state generation we analyze the evolution of the 
system, initially in the ground state, subject to a common 
bias pulse on the two gates, i.e., Vi = V2 = VoO{t) where 
0{t) is the unit step function. Then ([2]) reduces to 



H = hY^ cj[a]a^- + -(a^- + a])'^] + fia;(5(t)(a2 + 4) + ^i2, 

(4) 



where cj^ 



25i,2/M, Hi2 = FqUl e 

(t) 



Dh/M^cu^, and 6{t) = V^i^V^/VUM^. The af are 
defined through pj — ipo{cLj — aj)/V2, and qj = xo(a^ + 
aj ) / a/2 with xq = y/h/Muj and po = V Mfiuj. 

We consider the situation when the system is cooled to 
ksT <C fiuo and at time t = the flexural modes are in 



their ground states |0). Equal voltage pulses are applied 
to both gates 5{t) = So6{t)^ inducing system's evolution. 
As mode 1 will not be appreciably affected by the weak 
coupling term i^i2, the dynamics of the modes decouple. 
Hence, mode 1 will remain close to its ground state for 
t > 0. The remaining mode 2 describes a particle in a 
potential v{^) = \/2So^ + ^^/2 + e^^ with ^ = q2/xo and 
an equilibrium position ^o- Introducing the displaced 
oscillator operator 62 = a2 — ^o/>/2, and applying the 
rotating wave approximation (RWA), the Hamiltonian 
becomes H = TiCofih + 3hujehl/2^ where 



hb = blb2, 



(5) 



If e^o ^ situation is similar to the one in [9]. The 

intial state = 0)) = |0)^^ resembles a coherent state 
in the 62-basis, i.e, lO)^^ ~ |— ^o/V^)^^- If the system 
evolves a time Ti = 7r/(3cj£), we expect to find the state 

l^(Ti)) oc e--/4| - eoe-"'VV2>62 + e'^/^l^oe-'"' /V2)b, 

which, to an overall phase, is in a-basis given by 

IV(Ti)) cc |eo(l-e-^"*)/V2)a, +i|6(l + e-'-*)/\/2)a,. 

(6) 

To verify this we simulated the dynamics using the full 
two-mode Hamiltonian in ([2]) and a Hilbert space of 40^ 
number states in the occupation basis. The parameters 
used were L = 126 nm, uq = 60 nm, Tq = 0.003 N/m, 
and Vo = 0.23 V. This corresponds to a;/(27r) ^ 560 MHz, 
e = 7 ■ 10~^ and Sq = 1.3. The position shift becomes 
^0 - -1.8 ^ -V2S0. 

In Fig. ^ the probability density of finding mode 2 at 
a position ^ is shown as function of time. Only instants 
where t = 27Tn/uj are sampled hence fast oscillations with 
frequency uj are not visible. The initial state lO)^^ = 
I —^0 / V2)^ evolves into a first cat-like state at time Ti ^ 
0.42 /is ~ 7T/{3uje). This and the following cat-like state 
emerging at T2 ~ 1.27 ps ^ 7r/{uje) are marked with 
vertical dashed lines. The simulations also verified that 
mode 1 remains close to its ground state. 

To read out the state the two gates may form part 
of a capacitor in an LC-circuit with cjlc ^ co' in the 
resolved sideband limit. This allows for side band cooling 
and measurement of the quadrature operator [25[ Xjy = 
[e^^al + e~^^a2]/ V2. In Fig. [2]3 we show the envelopes 
of the average (Xj^) along with the quantum fluctuations 
(AX^) as functions of time. The fluctuations have the 
form (AX^) - A{t)^B{t) cos(2a;t), where A and B vary 
slowly in time. A signature of the built up cat state is 
the reduction of {Xjy) with an increase in (AX^). 

For the cat state to emerge, dissipation must be weak. 
Even at zero temperature damping will cause decoher- 
ence. As shown in [nl, [13 this occurs on a time scale 
tc ^ where Q is the resonator quality factor. Re- 

quiring tc ^ Ti yields Q ^ ^o/^- -^^^ protoype 
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FIG. 2: (Color online) (a) False color plots of snapshots of 
time evolution of the position probability distribution for one 
flexural mode. The snapshots are taken at the turning points 
of the corresponding classical trajectory of the system. The 
positions (y-axis) are scaled to the quantum zero point fluc- 
tuations xq. with vertical, white dashed lines, (b) Corre- 
sponding time evolution of the envelopes of the quadrature 
{Xu) (continous line) and its associated quantum fluctuations 
(AX^); A{t) (dashed line), B{t) (dashdotted line). The ap- 
pearance of a cat state is signalled by a decrease in (Xi,) 
along with an increased contribtion from quantum fluctua- 
tions (AX^) in the noise of (Xu). Vertical dashed lines indi- 
cate the agreement with Fig. 2a. 



system, this inequality is fulfilled if Q > 10^. Recently 
Q-factors up to Q ~ 10^ were reported in graphene res- 
onators [l6| . 

The cat state above is a product state between modes 
1 and 2 in a-basis. We now demonstrate a cat-like state 
involving both modes that is not a product state in this 
basis. Preparing the system in the ground state |0,0)^ 
and switching on only one gate at time t = 0, i.e., V2 = 0, 
Vi = the degeneracy is lifted. Just as the operators 



a^-^^ and b)^^ respectively diagonalize the linear part of 



when Vi,2 = and Vi = V2 7^ 0, the normal modes 
when ^1 7^ = are found by diagonalizing the linear 
part of (j2j) by means of the transformation 

di = cos{6)ai — sin(6>)a2, 

d2 = sin(6>)ai + cos(6>)a2 + r?o- 



Here 7^0 = k.V'^ /VhMuj^ and we have neglected terms 
of order {{Qi - (^2)/^^! < 1, where Qj^^ = cj^ T K12/M 



are the new eigenmode frequencies. The Hamiltonian (j2j) 
now transforms to 



(7) 



where Hnl is quartic in d-operators. As parameters like 
length and tension are never exactly equal in both x and 
^/-directions, the two modes 1,2 are only approximately 
degenerate, i.e. uoi 7^ co'2. This leads to the requirement 
1(1^1 - 9.2) /{oji - a;2)| > 1 so that 7r/4, and 



Hnl 



■d\?{d: 

f2, ,3 
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2%)^ 



(8) 
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e = h{2D - F) /AM^uj-^, A = h{QD - F) /2M^w^ 

The initial state 10,0)^^ is in cZ-basis |0,7?o)d- After evo- 
lution with (|7| this state will, analogous to the situa- 
tion studied above, enter a superposition at time Ti « 
7r/(3we) « 



l^(fi)) cc |0,%)rf + i|0,-r?o>rf + a|x). 



(9) 



Here a\x) = <^Xln=i \^)di I^^)d2 quartic 
coupling term in (|8]). The state component 1 0,770)^ + 
i |0, —rjo)^ corresponds to a cat-like non-product state 

|0,0)^ + z |-\/27/o, -V2r/o)^ in a-basis. 

To verify the creation of this state we numerically ana- 
lyze the evolution of |0, 0)a. Fig. 3a shows the Wigner dis- 

tributionipy(a,^i) = (2^)-^^^ ^^(6 +e/2|pi|6 - 
^/2) of the reduced density matrix pi of mode 1 in a- 
basis at Ti . The distribution is plotted as function of the 
dimensionless position ^1 and momentum tti and has a 
bimodal structure. The distribution for mode 2 is iden- 
tical, W{^2,t^2) = W{^i,t:i). 

To remove a\x) from ([9j) we introduce the projection 
operators = (|^) {'^\)d^ ^ study (Fig. 3b) the 

Wigner distribution of the projection Pq\iI){Ti)). One can 
here clearly recognize the distribution corresponding to 
the state |0, 0)^ + i | — >/27/o, — which is displayed 

in Fig. 3c for reference. To ensure that the a\x) compo- 
nent is not of major significance, the time evolution of 
the projections (Po(^)), and (/ - Po(^) - 

are shown in Fig. 3d. 

The expression for |V^(Ti)) in (|9|) is in RWA. In the 
Schrodinger picture the state has a fast oscillating com- 
ponent and is at t = Ti 



Vo 



|0, rjoe 



)d + i\0,-r]oe-'^'')d + a\x) 
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FIG. 3: (Color online) (a) False color plot of reduced Wigner 
distribution of the time evolved initial state |0,0)a, sampled 
at Ti as function of the dimensionless position and mo- 
mentum TTi. The reduced Wigner distributions are iden- 
tical for both modes VK(^i,7ri) = W {^2,712). A bimodal 
structure is seen. (b) False color plot of Wigner distri- 
bution of Poy9i(^i, TTi, Ti)Po, clcarly demonstrating bimodal- 
ity. (c) False color plot of Wigner distribution of |0,0)^ + 
i I — V^^o, — a/2^o)^- The similarity to the projection in (b) 
is evident, (d) Time evolution of projections (Po(t)), (P2(t)) 
and (/ — Po(t) — P2(t)). The most significant contribution 
comes from (Po(t)). 

Here Q2 is the eigenfrequency renormalized due to the 
nonlinear it ies in (j8]) [cf. (j5j)]. 

Similar behavior is seen for evolution from the initial 
state |?7o,0)^ = |0, — a/2?7o)^. Bimodality is then observ- 
able at T2 ~ 0.75/i5. We attribute the difference T2 — Ti 
to the position shift 2r]o in ([8]). 

Finally we demonstrate cat-like non-product states in 
both a— and tZ— bases. Assume the cat state (|6]) was gen- 
erated by the two-gate configuration at Ti. One gate is 
then switched off when CjTi — 27rm, m=integer. This 
kind of time-domain control has been shown possible 
in [Isl . The cat state is then a superposition of coherent 
states in d-basis |V^(Ti)) ex |0, |-^o, ^0 + ^o)^. If 

V = >/2Vb, then 7^0 ^ —^0, and the state is 

|^(ri))«|0,r?o), + i|%,0),. (10) 

As in previous cases one would expect that after an evo- 
lution with ([7]) both modes would enter a superposition 

|0, 7?o)^ + i |0, -7?o)rf + i (|7?0, 0)^ + i |-7?o, 0) J + /3|C), 

where p>\Q again is a small remainder due to the quartic 
coupling. Numerically we observe cat-like states in both 



modes in time interval 0.72/is < T\ < 0.82/is. These are 
cat-like superposition states which are not product states 
in either a- or d-basis. 

We have shown that due to the intrinsic non-linearities 
in graphene, generation of cat states and multimode cat 
states is possible by local back-gate manipulation. The 
nonlinearities are strong enough to avoid the necessity 
of coupling the membrane to an auxilliary system. To- 
gether with recent advancement in graphene device fab- 
rication and improvement of cooling schemes this opens 
for further fundamental studies of macroscopic quantum 
phenomena. 
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